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Methods are developed for calculating the normal coordinate vibrations of isolated helical homo- 
polymers and of the antiparallel sheet structures formed by some helical polymers in the solid state. 
The dynamical equations are expressed in Cartesian displacement coordinates starting from an internal 
coordinate harmonic force field. As an example of the method the dynamical equations of the honey- 
comb lattice are derived. 
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1. Introduction 

An essential feature of many polymeric crystals is 
layered structures which are held together primarily 
by interchain hydrogen bonding. Examples of poly- 
meric materials which form this type of structure 
include some nylons, synthetic polypeptides, and 
fibrous proteins. Interchain bonding is expected to 
influence vibrational properties especially in the low 
frequency region and thereby affect the infrared and 
Raman spectra, the inelastic neutron scattering, and 
thermodynamic properties such as the low temperature 
specific heat. 

With the availability of commercial Fourier trans- 
form infrared spectrophotometers and low light 
scattering Raman spectrometers low frequency vibra- 
tional studies on polymeric systems are now practical 
and several such studies have recently appeared in 
the literature [1-3]. The assignment of vibrational 
spectra obtained from solid state samples has been 
based mainly on normal coordinate calculations of 
isolated chain models [4-7]. Interchain effects have 
been either totally ignored or treated in part by pertur- 
bation methods [4, 8, 9] because the vibrational secular 
equation may become quite large. For example, the 
normal mode frequencies of polyglycine, (CH 2 NHCO) w , 
in the fully extended conformation, polyglycine I, 
may by calculated from an isolated chain model with 
1 C 2v symmetry by solving a secular equation of dimen- 
sionality twenty-one. For the more realistic two- 
kdimensional pleated sheet model the dimensionality 
Eof the secular equation is 84 since four glycine residues 
"comprise the two-dimensional unit cell. To calculate 
/ the vibrational frequency as a function of the phase 
differences (phonon dispersion curves) requires the 
I solution of a secular equation which is twice as large 
because the matrix elements are, in general, complex. 



The phonon dispersion curves are useful for com- 
parisons with the inelastic neutron scattering data 
[7, 10] and with the low temperature specific heat [7]. 
In this paper a method is introduced for calculating 
the normal coordinates in terms of the Cartesian dis- 
placement coordinates of isolated helical homo- 
polymers. This method is extended to antiparallel 
sheet structures with four chemical repeat units com- 
prising the two-dimensional unit cell. The method is 
demonstrated for the honeycomb lattice shown in 
figure 1. 

2. Method for Calculating the Normal 
Coordinates of an Isolated Chain 

The adaption of the Wilson GF method [11] to the cal- 
culation of the normal coordinates of helical infinite 
polymers was first described by Higgs [12] who utilized 
symmetry adapted internal coordinates. Miyazawa and 
coworkers [13] extended Higgs' method by defining the 
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Figure 1. Honeycomb lattice unit cell is indicated by dotted lines. 

The polymer chains lie along the z axis. The interchain bonds are parallel to the x axis. 
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internal coordinates in terms of the helical parameters. 
This method is rather cumbersome especially for com- 
plicated chemical repeat units. Small, et al. [6] have 
developed a method by which the vibrational secular 
equation may be set up from the internal coordinates of 
one chemical repeat unit defined in the Cartesian co- 
ordinate system. The dynamical matrix in terms of 
the internal coordinates in nonhermitian so as to re- 
quire two separate matrix diagonalizations. Further- 
more, there are usually redundant internal coordinates 
introduced so that the dimensionality of the secular 
matrix is larger than 3M where M is the number of 
atoms in the chemical repeat unit. Both of these dis- 
advantages are removed by expressing the vibrational 
problem in Cartesian displacement coordinates. An 
expression for the vibrational secular equation for 
helical polymers in the symmetrized Cartesian dis- 
placement coordinate basis set will now be developed 
following the method devised by Small, et al. [6], for 
the internal coordinates. 

The kinetic energy in terms of the Cartesian dis- 
placements is 



TV/2 



2T= y y X r mX r 
^ Zj j j j 

r=-A72+l j=l 



(1) 



where Xr is the jth displacement coordinate of the 
rth chemical repeat unit for a polymer of N total repeat 
units with reentrant boundary conditions and rrij is 
the mass of the atom undergoing theZr displacement. 
The potential energy assuming harmonic forces is 



where 



2V=^^Xrfr ; PXf 
r, P j, I 

fj>p=d 2 VldXjdXf 



(2) 



are force constants in terms of Cartesian displacement 
coordinates. 

In matrix notation eqs (1) and (2) become 



and 



2T=^X r MX r 

r 

2V=^ P£' p lk p 



where the tilde means the transpose vector, 
potential energy, eq (4a), may be rewritten as 

2V=y V Xrfr,r+sXr+s 



(3) 
(4a) 
The 

(4b) 



where the index s labels the neighboring units of the 
rth unit and£ r + s = for \s\ > k. 

The symmetry coordinates, X(0), are linear com- 
binations of equivalent coordinates. Unlike the internal 
coordinates which are invariant of unit location and 
equivalent, the Cartesian coordinates depend on unit 



H r (i//) is comprised of M 3X3 blocks along the 
Hlagonal where each block is of the form 



(cos(n 
sin (n 
> 



(n/0 
(nil) 



-sin(n//) 

cos(nfj) 




?) 



(5) 



and i// is the helical angle. 

The matrix H r (i//) rotates the x,y coordinates of 
the Oth unit into those of the rth unit. The form of trans- 
formation given in eq (5) assumes a right-handed helix 
with helical axis Z and with r > labeling units reached 
by clockwise rotations plus translations. 

The normalized symmetry coordinates X(0) arc 
defined as 

X(0) = l/VN J) exp(-iV0)H- r (</>)X>\ (6) 

r 

Solving eq (6) for X r and substituting into eq (4b) yields 
2^=2 2 exp[ir(0' - 6)] 

r,s 8,6'. 

X(0)U- r Wi r ' r+s U r WU s WMO f )exp (isO') (7) 

From the definition of the transformation matrh 
H r (i//) and the covariance of £ r > r + s it follows that 

H- r (</>)£ r > r+s H r (i//) =£ ) ' s . 

With this result eq (7) becomes 

2V= 2 X(0)i^E s WX(0) exp (isO). (8) 
s,e 

The potential energy may also be expressed in the 
internal coordinate basis set as 



2V=^WF^R9 

t,q 



(9) 



where R 9 is the column vector of internal coordinates 
in the q\h chemical repeat unit and 

F\f= dWIdR^R* 

are the internal coordinate force constants. 

The internal coordinates are defined through tl 
B matrices [6], 

p 
so that eq (9) becomes 

2F=2 2 % r B t$r E ttq B z qtP X p 

t, q p,r 

which yields upon comparison with eq (4a) 



location 
through 


but 
the 


may be made to form an equivalent set 
rotational operator jj r (i//). The matrix 


fr,p = 


= 28' 

t,q 


,r^_t 


,qftq,p. 1 
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Utilizing the expression [6] 

and the fact that E/»« = E ^-' we have 

£0,* = ^ H^i//)^ '-^ ^ ^-^!^)^^^ (10) 

t, r 

Where the summations are over positive and negative 
values. We find by substitution of eq (10) into eq (8) 
the potential energy in terms of the Cartesian displace- 
ment symmetry coordinates and the internal coordinate 
force matrix: 

2F= S E X(^)H-H^)B°^exp(-^)F ^ 

exp (i>0)B°'*B^)exp(;s0)X(0). (11) 



Defining 



and 



E(0) = X exp U's0)B°' s H*(i//) 
F(0) = ^ exp(i'r0)F°-'\ 



(12) 



(13) 



Equation (11) becomes 

i 2F=^X(0)B(0)F(0)B(0)X(0). (14) 

e 

The kinetic energy, eq (3), becomes in the Cartesian 
symmetry coordinates 

2T=^X(0)MX(0) 
e 

or in the conjugate momenta 

2T=^P(d)M- l P(0). (15) 



From eqs (14 and 15), we find through Hamilton's 
equations the eigenvalue equation 

M- 1 B(e)E(fl)B(fl)L^(e) = Lr(fl)A(9) (16) 

here A(0) is the diagonal matrix of the squares of 

^e eigenfrequencies and L x is the matrix of eigen- 

ectors. Equation (16) may be symmetrized in the 

irm 

i- 1/2 B((9)E(0)B((9)M- 1 / 2 L;(0)=L;((9)4(0) (17) 

/here L^(0) —M 1,2 Lix is the matrix of eigenvectors 

n the mass-reduced symmetry coordinates. 

The form of eq (17) is identical to that given pre- 

usly [14] for translational symmetry (i.e. H r (i//) 

the unit matrix). The essential feature in this deri- 

tion is that the dimensionality of £(0) equals the 

mber of independent coordinates in the chemical 
peat unit and that only the internal coordinates 



of one chemical repeat unit in terms of Cartesian dis- 
placement coordinates need be specified. 

3. Two-Dimensional Lattice of 
Antiparallel Chains 

Consider a two dimensional lattice, figure 1, in the 
X, Z plane with two chains passing through the unit 
cell and helical axes parallel to Z. The helical angle, 
i//, is 180° and there is a two-fold screw rotation axis 
along X. The latter symmetry operation reverses the 
sense of the two chains to produce the antiparallel 
structure. The translational unit cell contains four 
chemical repeat units numbered (0, 0), (0, 1) in the 
a. chain and (1, 0), (1, 1) in the /3 chain. 

The four units are related by the two-fold screw 
axes along Z and X. For the two-dimensional lattice 
we have [14] 



^here 



f(4>, 0) =§(<£, 0)F(</>, 0)B(c/>, 0) (18) 

B((/>, 0) = ]? W s ex P [i(rcf> + sV) \ 

r, .v 

F((/>, 0) = 2 F r ' s exp [i(n/> + s0)] 



with 



= k z -tz. 

The quantities tx, tz are the lattice vectors directed 
along the X and Z axes and k*, kz are the wavevectors. 
The quantities (f) and are the phase differences 
between identical displacements in neighboring unit 
cells. The dimensionality of the matrix £((/>, 6) is 
3/n X 3m where m is the number of atoms in the two- 
dimensional unit cell whereas the dimensionality of 
F(</>, 0) equals the number of internal coordinates 
which is usually greater than 3m since redundant 
coordinates are included. 

The basis set used to express the quantities in 
eq (18) are linear combinations of equivalent chemical 
repeat units: 

, N x-i N z-i 

Xj,i(<l>,0)=llVNxTz 2 2 



r = * = 



exp [— i(r<j>- 



td)}X r > 



h i 



(19) 



where j, 1 = 0, 1 label the chemical repeat unit and 
r, t the unit cell. The symmetry coordinates, eq (19), 
may be related to those of eq (6) as follows. The 
chemical repeat units are renumbered so that X 7 '' q 
equals Xjfr j)/2 ' (9_/)/2 with 7, / = or 1 forp, q even or 
odd. A generalization of eq (6) yields the symmetry 
coordinates 

2Nx - 1 2Nz - 1 

x(<t>',e') = ii2VN x N z jr £ 

l> = q = 

exp [-i{p(f>' + qe')]H-i>' ~ "X"- « (20) 



353 



464-124 OL -72 - 3 



where H p ' 9 is the transformation matrix representing p 2-fold rotations about X and q 2-fold rotations about 
Z. The matrix jj y ' - q is comprised of m 3 X 3 blocks with each block given by 



-1)" 

(-1)* + * 





S). 



Using the relationship between Xf>Q an d Xj'p eq (20) may be rewritten as 
X (0' , ft') = 1/2 VA^VT ( 2 2 exp [- £ (p0' + q6')] Xg§ " /2 

l/>, even q, even 

+ 22 exp [-i(p</)' + ^')]H»>'Xg/v ( «- 1)/2 

l>, even g, odd 

+ 2 2 exp f -8 ' (P*! + «»')] Hi.«XiPo 1)/2 ' 9/2 

/>, odd <?, even 

+ 2 2 ex P t~* (P*' + fl*')] g^X^/*. <€-»>/« 

Pi odd 9, odd 



(21) 



X (</>', 0') = 1/2 [X 0> „ (20',20') + H o,, X D> , (20',20')exp(-tft') 

+ H 1 ' < Xi,o (20', 2ft') exp (-i«p') + H 1 - 'X,, , (20\ 20') exp(-i(ft' + 0'))]. 

Functions orthogonal to X(0. 6*) are X(0. ft+ 77) . X(0 + 77. ft) . and X(0 + 77. 0+ 77) which are related to those 
ofeq(19) 



X (0, 0) 




I 


e -,»H«. 1 


e-WRi.o 


e -i(*+9)JJl, 1 




Xo,o (20,20) 


X(0, ft+77) 


1 


I 


-e-'OH - 1 


e-^H 1 - 


_ e -/(«+fl)JJl, 1 




Xo.i (20,2ft) 


X(0 + 77, ft) 


2 


I 


e-' e H°- » 


_ e -i*gi,o 


_ e -i(rf>+6l)JJl, 1 




Xi.o (20,20) 


_X(0+77, ft+77) _ 




I 


-e-^H - 1 


-e-^H 1 - 


e -««+8)gl, 1 




_X,, 1 (20,20)_ 



(22) 



or 



X (0, (9) 

X (<M+7T) 

X (0+77, 0) 

X (0 + 7r,0 + 7r) 



u(<M) 



X o ,o(20, 20) 
Xo,i (20, 20) 

Xi,o (20,20) 
Xi,i (20,20) 



(23) 



where I is the identity matrix. 

We now derive an expression similar to eq (22) for the internal coordinates so that £(0, 0). eq (18). may 
be transformed to the new symmetry coordinates. It is convenient to list the internal symmetry coordinates 
of the lattice as [R ,o (0. 0), So, i (0. 0), r a (0. 0). R,.<> (0. 0). Ri.i (0. 0), r^ (0, 0)] where R } , (0,0) 
are intrachain coordinates of the 7, / repeat unit. r a (0. 0) and r $ (0. 0) are interchain coordinates associated 
with the a and ft chain respectively. In a manner analogous to that which lead to eq (22) we find 
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R(<M) 

R(</>, 0+tt) 

r(</>, 9) 



R(0 + tt, 6/) 

R((/>+7T, 0+77") 

r(0 + 7T, 0) 



. c w 








-i<£ 



-e !W 




V^ 



-e- 

-e- 

V2 



<</> 



;</> 



e -!'(<* + «) 





_ e -i(*+e) 
e - «* + «) 









V2« 








- V2 e -«• 



R o> o(20, 20) 
Ro,i(20,26) 
r„(20, 0) 



R,,,(20, 20) 

R,,,(20, 20) 
rp(20 0) 



(24) 



W(0, 0) 



(25) 



R„,o(20, 20) 
Ro, 1 (20, 26/) 
r a (20, 0) 

R,,o(20 20) 
«i.i (20,20) 
^(20, 0) 



The H'''' matrices do not enter into W(0. 0) since the internal coordinates are invariant to position within 
the unit cell. Using eqs (23). (25). f (0. 0) is transformed to 



where 
and 



f'(0, 0)=B'(0, 0)F'(0, 0)R'(0,0) 
B'(0, 0)= W(0, 0)B(20, 20)0(0, 0) 
F'(0, 0) = 1(0, 0)F(20 20)0(0, 0). 



(26) 



In the absence of interchain interactions £' (0. 0) consists of four noninteracting blocks, two of which 
correspond to intrachain phase difference and the other two to 0+7r. It will be shown that with certain types 
of interchain interactions £ (0, 0) contains nonzero elements only between the blocks (0, 0) and (0 + 77,0+77) 
and between (0, 0+ 77), (0 + 77, 0) so that the dimensionality of the dynamical matrix has been reduced by a 
factor of two. 

Let 



b; (0, 0) = 1/V2 



1 



1 



-./. 



B«(20,20) 



Bf (20, 20) 



0(0,0) 



(27) 



where B? (20, 20), B? (2(/>. 26) are the B matrices associated with coordinates r a , rp respectively. These B 
matrices may be expanded in terms of contributions from each of the four residues in the unit cell. 
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B«(2</>, 20) 




B a 

=0,0 


B a 

=0, 1 


&?,« 


I"! 


Bf(2*,2tf) 




B^ 

=0,0 


B^ 

=0, 1 


21. 


Bf fI 



Letting, 



we find upon substitution of eq (28) 



B«(0, fl)=B« (20,20) U (0,0) 



1 



and 



B«(</>, 0)=- S exp[i(j(t> + kd)]Bj k Yi-J>- k 
2j, /t = o 



B*(<M)=^ V exp[*O*<£ + *0)]B* H-i-*. 



Substitution of the expression for B"(20, 20)lJ((/>, 0) into eq (27) yields 



b; (*, e) = 



where 
and 



B+((/>, 0) B + (0, 0) B + (0, 0) B + (^, 0) 
B -(0, 0) B-((/>,0) B"(0, 0) B-(0,0) 
B± (0, 0) = 1/ V2 [B a (0, 0) ± exp (- j0)B* (0, 0) ] 



0=0+ 7T, </>=(/> + 77. 

If a diagonal valence force field is assumed for the interchain interactions 



F'(<M) = 



F(0) 














°\ 





FW 











° \ 








F° 




















F(0) 




















F(0) 


0/ 

















F / 



(28) 



(29a) 



(29b) 



(30) 
(31) 



(32) 



where F(0) is the isolated chain force matrix, eq (13), and F is the valence field force matrix for interchain 
interactions. From eqs (30-32) we have 



f ' (4, B) 



MB)+M<t>+6) /.(<M;<M) M<t>,0;$,e) M4>,e-,$,9) 

f r (6)+M<f>,0) M<f>,O~;$,0) f,(4>,d;$,V) 

cc. MB) +fi($, 9) M$, 0; $, 6) 

f r (S) +fi($, 3) 



(33) 



where f r (6) is the symmetrized force matrix for the isolated chain [from eq (16), f r (B) = B(0)F(0)B(0)], 
M4>, B; <*>, B) = B + (</>, B) F°B + (</>, B) + B (<j>, B) F°B (</>, B) (33b) 
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and 

The vibrational secular equation for the two-dimensional lattice is 

| M-^f' (0, ^M-^-co 2 ! | = 0. 



(34) 



A relationship exists between B a ((/>, 0) and B^ (0, 0) which is helpful in the further simplification of/' (</>, 6). 
The quantities B a A . ($, 0), B^ (0, 0) which appear in eqs (29 a. b) are sums of B matrices multiplied by phase 
factors. 

fi«, (20, 20) = V exp [2i(l+ + nO) £«;/•" (35a) 

and 

Bf >Jt (20, 20) = ^ exp [2i (Z0+ "0)]Bf;^ (35b) 

/, n 

where the superscripts /, ra label the unit cell. As the internal coordinates given byB? fc (</), 0) and B" ; . (</>, 0) 
are equivalent it follows that 

and from eq (35b) 



j,* V 



=j. fc 



R/3,/,« = R a ' '+i-*» -" H 10 



8,% (20, 29)= exp [~2i (/-l)*] &%,!-* (2*, -20) H ' ° 
Substituting the above expression into eq (29b) yields 

jf(</>, 0) = exp [»(0+ 0)]B«(0, - 0)H°- » 
and from eqs (31) and (33b) 

/(<£, 0;<£\0')=B a ((/), 0) F°B«(c/>\0')+exp [-* ((/>-</>' + 0-0')] H ' -* 

B« (</>,- 0) F°B« (</>', -0') H 01 . (36) 

The secular equation may be constructed from the B and F matrices for the isolated chain and the inter- 
actions between the two chains in the unit cell. Under certain conditions as shall now be shown f ' (</>, 0) reduces 
to block form. 

We assume the polymer possesses a plane of symmetry and that the two chains in the unit cell are coplanar. 
If the internal coordinate connecting the two chains is a bond stretch directed along the x axis as in figure 1 
B a (<£>, — 0) is independent of and 



fi (</>, 0;4>\0') = [14- exp [-i (</>-</>' + 0-0')]] B« (</>, 0) F°B« (<£\ 0'). 



(37) 



The quantity £/ (<f>, 0;<£\ 0') =0 if(<f>', 0') = (<£ + tt, 0) or ($', 0') = (<)> 9 + ir) and T(0, 0), eq (33a) becomes 
upon rearrangement 



f'(<M) = 



fr(o)+fi(<t>,e) 


fi (0,0; <M) 










f r (lf)+Mj,V) 








c.c. 




/r(9)+/,(<M) 


/K<M;$,0) 








/r(«)+/((?,«) 



(38) 
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Hence, for interchain interactions directed along the x axis (perpendicular to the chain axis) the chain 
modes given by 0, interact with those of the other chain in the unit cell given by c/>-f7T, + 7T with the inter- 
action given by 2g «((/>.#) E () Br(</> + 7r, # + 77 ")- 

A dynamical matrix of the form of eq (38) is retained for interchain interactions which couple x and z mo- 
tions if the chain possesses local C 2 symmetry about the interchain bond as in the case of the honeycomb lattice, 
figure 1. Consider the internal coordinates which are deformations of the angles formed by the interchain bond 
in figure 1. The two angles with common apex atom are equivalent because of the local Co symmetry. Under the 
two-fold screw rotation about x these internal coordinates are carried into an equivalent pair in the second 
chain. These four internal coordinates are therefore equivalent and if the equivalent ones in I|f are renumbered 
it follows that 

, n+k 



from which it can be shown that 






B# (</>, 0) = exp [i ($-0)] B« (0, 0) H 1 



For in-plane internal coordinates 



and eq (38) follows. 



B^ (0, 0) =- exp [i (0-0)] B a (</>, 0) 



4. Example: Honeycomb Lattice 



The preceding method will now be applied to the honeycomb lattice shown in figure 1. The in-plane intra- 
chain internal coordinates of the honeycomb lattice are /?" bond stretch and /?". angle bend; one of the set 
of two equivalent interchain coordinates consists of r". bond stretch and r 2-5 . angle deformations. Let 2e be 
the intrachain angle. 5= sin e, c=cos e and r, R the intrachain and interchain bond distances respectively. 
The intrachain coordinates are 



\R$J \slr-c/r/— \-2slr0 J— \s/rclrj — 



and the interchain coordinates 











0\ /-I 

\ I sir -(clr+HR) 

X?;r'+ sir cJr+UR 
sir c/r I HR 

0/ \0 -HR 



Xo, o 4 
_o, l ^ 






o\ 








sir 


-c/r 











0/ 




XO, i 
_1,0 ' 



/o 







is/r 



a \ 







-c/r 



XO, 
1, 1 



From eqs (12 and 13) the dynamical matrix for the isolated chain is. assuming a valence force field. 

cos 2 0/2(/ r c 2 + 4 5 2 / € /r 2 cos 2 0/2) 2isc sin (9(/ r -4/ e /r 
Use sin Oifr-tyjr 2 cos 2 (9/2) 4/ r 5 2 sin 2 0/2 + 4/ e , 
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[4 cos 
r 2is 



fjr 2 cos 2 0/2)1 
if € lr 2 c 2 sin 2 0J 



where f r , ft are the bond stretch and angle deformation force constants respectively. The interaction matrices 
are found through eqs (28). (29a). and (37) to be 



/((M;<M) = 



yt<M;<M = 



"2/fl sin 2 (0/2) +8/«5 2 /r 2 cos 2 0/2 -4*/ 5 s/r sin (6) [c/r+ cos 2 (<f>/2)/R] 

4i/ 8 s/rsin (tf)[c/r+ cos 2 (0/2)//?] 8/ 8 (c 2 /r 2 sin 2 ((9/2) + l//? 2 cos 2 (4>/2)+2c/rR 

cos 2 (</>/2)sin 2 (<9/2)) 



i/ft sin 
_—2f$s/rR sin sin 



■2fhSJrR sin </> sin 
-4i/ 8 //? sin (*)(1/J? + c/r) 



where ./«,/« are the interchain force constants for bond stretch and angle bend respectively. From /(</>. 0; 0. 0) 
and /(0. 0; 0. 0) the first block of f '((/>. 0). eq (38), may be constructed. The resultant secular equation may 
be solved for values of ^ </> ^ 7r/2 and ^ 6 ^ 77 to obtain the phonon dispersion curves over the z-extended 
zone. The variation of phonon energies with k c for a particular value of hj are found by folding the phonon 
dispersion curves about = 90°. 



5. Summary 

Utilizing the unit cell symmetry, the dynamical 
matrix (eq (33a)) of the antiparallel chain lattice has 
been derived. The elements of this matrix are ex- 
pressed in terms of the dynamical matrix for the 
isolated chain and the B and F matrices for the inter- 
chain coordinates. The dynamical matrix of the isolated 
chain and hence that of the lattice may be determined 
by specifying the internal coordinates of one chemical 
repeat unit. It was found that with certain types of 
interchain forces the dynamical matrix of the anti- 
parallel chain lattice may be reduced to block form 
for arbitrary values of the intrachain and interchain 
phases. All normal modes may be determined by 
solving one of the two blocks of resulting secular 
equation for ^ <\> ^ 7r/2 and 0^6 ^tt. The reduction 
in the dimensionality of the secular equation with 
certain types of interchain forces can lead to a sig- 
nificant reduction in computation time especially for 
more complicated homopolymers such as polyglycine. 

The methods used to derive the dynamical matrix 
of the antiparallel chain lattice may also be applied 
to other polymeric crystal structures. 
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